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ABSTRACT 

The physical conditions within damped Lya systems (DLAs) can reveal the star for- 
mation history, determine the chemical composition of the associated ISM, and hence 
document the first steps in the formation of present day galaxies. Here we present cal- 
culations that self-consistently determine the gas ionization, level populations (atomic 
. fine-structure levels and rotational levels of H2), grain physics, and chemistry. We 

' show that for a low-density gas {tih^ 0.1 cm~^) the meta-galactic UV background 

If^ , due to quasars is sufficient to maintain H2 column densities below the detection limit 

■ (i.e A^(H2)^ 10^** cm"^) irrespective of the metallicity and dust content in the gas. 

Such a gas will have a 21 cm spin temperature in excess of 7000 K and very low C i 
and C II* column densities for H i column densities typically observed 50 per cent in 
DLAs. 

O I We show that the observed properties of the ~ 15 per cent of the DLAs that do 

show detectable H2 absorption cannot be reproduced with only the quasar dominated 
^ ' meta-galactic UV radiation field. Gas with higher densities {tih^ 10 cm~'^), a moder- 

ate radiation field (flux density one to ten times that of the background radiation of 
the Galactic ISM), the observed range of metallicity and dust-to-gas ratio reproduce 
all the observed properties of the DLAs that show H2 absorption lines. This favors the 
rS , presence of ongoing star formation in DLAs with II2. 

■ The absence of detectable II2 and C i absorption in a large fraction of DLAs can 

be explained if they originate either in a low-density gas or in a high-density gas with 
a large ambient radiation field. The absence of 21 cm absorption and C 11* absorption 
will be consistent with the first possibility. The presence of 21 cm absorption and 
strong C II* without II2 and C i absorption will suggest the second alternative. The 
iV( Al 11) / A''( Al III) ratio can be used to understand the physical properties when only 
C II* absorption is present. We find nn in components that show C 11* (without II2) 
is less than that typically inferred from the components with II2 absorption. We also 
calculate the column density of various atoms in the excited fine-structure levels. The 
expected column densities of O i* , O i** , and Si 11* in a high-density cold gas is in the 
range of 10"'^^ — 10"'^^ cm~^ for log N(H i)^ 20 and the observed range of metallicities. It 
will be possible to confirm whether DLAs that do not show II2 originate predominantly 
in a high-density gas by detecting these lines in very high S/N ratio spectra. 

Key words: Quasars: absorption lincs-ISM: molccules-Intergalactic medium- 
Galaxies: formation 



1 INTRODUCTION 

Damped Lya systems seen in QSO spectra are characterized 
by very large H I column densities, A'(H l)> 10^" cm~^, that 
are similar to those observed through gas-rich spiral galax- 



ies. The importance of DLAs in the paradigm of hierarchical 
structure formation can be assessed from the fact that the 
mass density of baryonic matter in DLAs at Zabs"^ 3 is sim- 
ilar to that of stars at present epochs (Wolfe 1995). Studies 
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of Lya and UV continuum emission from galaxies associated 
with DLAs usually reveal star formation rates(SFR) (or up- 
per limits) of a few Mq yr"^ (Fynbo et al. 1999; Bunker et 
al. 1999; Kulkarni ot al. 2001; MoUcr et al. 2002 and 2004; 
Weathcrloy ct al. 2005). The metallicity and depletion fac- 
tor in DLAs arc usually estimated from N{Zn ll)/Af(H l) 
and N{Fe ii)/iV(Zn ii) respectively (Lu et al. 1996; Pettini 
et al. 1997; Prochaska & Wolfe 2002; Ledoux et al. 2002a and 
Kharc ct al 2005). The inferred motallicities typically vary 
between log Z = —2.5 to for 2 ^ Zabs^ 3 with a median 
of ~ —1.3 (Ledoux et al. 2003). The measured depletions 
range between to —1.6 with a median value of —0.3. If 
dust content is defined as k = 10[z°/h] (i _ ioP«/z°l), 
then the median dust content in a typical DLA is k — 0.07. 
This is less than 10 per cent of what is seen in the Galac- 
tic ISM for a similar neutral hydrogen column density. If 
the physical conditions in DLAs are similar to those in our 
Galaxy or the Magellanic clouds then H2 molecules should 
be detectable. There have been very few detections of H2 
in DLAs (with 3 x lO" < iV(H2)(cm-2) ^ 3 x 10^®) despite 
extensive searches (Ge & Bechtold 1999; Srianand, Petit- 
jean & Ledoux 2000; Pctitjcan, Srianand & Ledoux 2000; 
Ledoux, Srianand & Pctitjcan 2002b; Levsliakov et al. 2002; 
Ledoux, Petitjean, & Srianand 2003; Reimers et al. 2003). 
Roughly 80 per cent of DLAs do not have detectable H2 
(with iV(H2) ^ 10" cm-2). 

The physical conditions in the H I gas can be probed 
by using the fine-structure absorption lines produced by the 
excited atomic species and the 21 cm transition. Apart from 
a few rare cases (for example Srianand & Petitjean 2001), 
C I is detected only in systems that show H2. The derived 
total hydrogen density (nn) based on the fine-structure level 
populations of the heavy elements are usually large (^20 
cm"^) (Ledoux et al. 2002b; Srianand et al. 2005; Wolfe et 
al. 2004). Like C I, C n* absorption is detected in every case 
where one detects H2. However, it has also been seen with 
lower column densities in a considerable fraction of DLAs 
without H2 (Wolfe et al. 2003a; Srianand et al. 2005, Wolfe 
et al. 2004). 

The search for 21 cm absorption in DLAs at Zabs^ 2 
have mostly resulted in null detections with typical spin tem- 
peratures ^ 10-^ K (Table 3 of Kanckar & Chengalur (2003) 
and Table 1 of Curran et al. (2005) for a summary of all 
the available observations). There are 8 DLAs at Zabs^ 1.9 
for which redshifted 21 cm observations are available. Red- 
shifted 21 cm absorption is detected for only two systems 
{zabs= 1-944 toward PKS 1157-^014 (Wolfe et al. 1981) and 
Zabs= 2.0394 toward PKS 0458-02 (Briggs et al. 1989)). The 
measured spin temperatures are 865±190 K and 384±100 K. 
However, none of these systems show detectable H2 (Ledoux 
et al. 2003; Ge & Bechtold 1997). The DLA toward PKS 
1157-1-014 is special as the QSO shows broad absorption 
lines and the Zabs of the DLA is close to the Zem of the 
QSO. The physical conditions in this system may not be a 
good representative of the general DLA populations. Inter- 
estingly, H2 is seen at Zabs= 2.8110 toward PKS 0528-2505, 
while no 21 cm absorption is detected in this system (Carilli 
et al. 1996; Srianand & Petitjean 1998). The lower limit on 
the spin temperature derived for this system is 710 K. How- 
ever, the excitation temperature derived from H2 rotational 
levels is < 200 K (Srianand & Petitjean 1998; Srianand et 
al. 2005). This system is also special since Zabs>Zem ■ The 



radiation field of the QSO has a much stronger influence on 
the physical conditions in this DLA (Srianand & Petitjean 
1998). The upper limits on the spin temperature derived 
for the rest of the systems are higher than 1000 K. The H2 
content of these systems is not known. 

Even though various properties of DLAs, listed above, 
have been investigated in detail (Matteucci et al. 1997; 
Prantzos & Boissier 2000; Howk & Sembach 1999; Izotov 
ct al. 2001; Vladilo, 2001; Liszt 2002; Hirashita et al. 2003; 
Wolfe et al. 2003a,b; Calura et al. 2003; Wolfe et al. 2004), 
very few attempts have been made to investigate all of them 
in a single unified calculation. That is the main motivation 
of this paper. 



2 CALCULATIONS 

The main aim of our study is to investigate the physical 
conditions in high-2 DLAs. In particular, our goal is to un- 
derstand the equilibrium abundance of H2, the excitations 
of H2 and atomic fine-structure levels, the ionization state of 
the gas, and the 21 cm optical depth, sclf-consistcntly. In the 
Galactic interstellar medium (ISM) H2 is usually detected 
either in a diffuse medium irradiated by the UV background 
radiation field or in high-density photo-dissociation regions 
(PDRs) near OB stars. One can anticipate this to also be the 
case in DLAs. At high redshift the diffuse UV background 
from QSOs will be an additional source of UV radiation. 

Recently there have been three attempts to model H2 in 
DLAs. Liszt (2002) uses two phase models similar to Wolfire 
et al. (1995) for this purpose. The models consider dust-free 
gas, so only the slow gas-phase H~ -I- H ^ H2 -I- e forma- 
tion process is important. The second attempt by Hirashita 
et al. (2003) models DLAs as large protogalactic disks. The 
density and temperature in the gas are determined by "non- 
linear hydrodynamic" effects. The radiation field is assumed 
to have a negligible contribution to the temperature of the 
gas and is used only for destroying the H2 molecules. Their 
model provides insight into the spatial distribution of H2. 
The third attempt by Hirashita & Fcrrara (2005) uses sim- 
ple H2 equilibrium formation models to study H2 in DLAs. 
Unlike Liszt (2002), no attempt is made to model ionization 
conditions of the gas and excitations of the fine-structure 
lines in the later two studies. The main aim of our paper 
is to use full self-consistent numerical calculations to under- 
stand (i) physical conditions in DLAs with H2 detections, 

(ii) the reason for the lack of H2 in most of the DLAs (iii) 
the origin of C II* absorption frequently seen in DLAs and 

(iii) the absence of detectable 21 cm absorption at high red- 
shifts (i.e z^ 2). 

The availability of good quality observational data al- 
lows us to estimate the metallicity, depletion, H2 abundance, 
H2 excitation, and populations of fine-structure levels in C I 
and C II (Ledoux et al. 2003; Srianand et al. 2005). One can 
hope to build more realistic models with all these in hand. 
This forms the prime motivation of this work. We study the 
ionization state, chemical history, and temperature of the 
gas using version 96 of Cloudy, last described by Ferland 
et al.(1998), and available at http://www.nublado.org. 
The details of the new H2 model are given in Shaw et al. 
(2005; hereafter S05). A comparison between predictions of 
our code and several independent calculations of PDRs can 
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be found at http://hera.phl.uni-koeln.de/~roellig. A 

direct application of a PDR is given in Abel et al. (2004; 
hereafter A04). The calculations presented here take into 

account various heating (e.g. dust photo-electric heating cos- 
mic ray heating etc) and cooling processes (for details see 
Ferland et al. 1998 and 805). 



2.1 The micro-physics of grains : 

We use the improved grain physics and molecular network 
as described in van Hoof et al. (2004), Ferland et al. (1994; 
2002), and A04. The grain size distribution is resolved, and 
the formalism described by Weingartner & Draine (2001a; 
van Hoof et al. 2004) is used to determine the grain charge, 
temperature, photoelectric heating, and drift velocity self- 
consistently with the local radiation field and gas tempera- 
ture. The extinction curves of grains in DLAs are not well 
known. We use a grain size distribution that fits the extinc- 
tion curve of the difi^use interstellar medium with Ry = 3.1 
(Table 1 of Weingartner & Draine (2001b)). We emphasize 
that the physical treatment of grains, and their effects on the 
surrounding gas, is fully self-consistent, and does not rely on 
general fitting formulae, such as those in Bakes & Tielens 
(1994). The grain charge and temperature are determined 
by the local gas conditions (mainly the electron density and 
temperature) and radiation field (including the attenuated 
incident continuum and emission by the surrounding gas, 
mainly Lya ). The result is a grain temperature and charge 
that depends on grain size and material. The temperature 
then determines the rate of H2 formation on grain surfaces - 
we adopt the temperature-material-dependent rates given 
by Cazaux & Tielens (2002). The charge establishes the 
floating potential of the grain, which then sets the grain 
photoelectric heating rate. 



2.2 Molecular hydrogen : 

The detailed treatment of the micro-physics of H2 is de- 
scribed in 805. Here we will briefly mention some of those 

processes. 

Generally, H2 forms via grain catalysis in a dusty cold- 
neutral gas. We use the total formation rate given by Cazaux 
& Tielens (2002) along with the size and temperature re- 
solved grain distribution described in van Hoof et al.(2004). 
This is an exothermic process and H2 is formed in exited vi- 
brotational levels, a process referred to as formation pump- 
ing. The results presented below use the state-speciflc for- 
mation distribution function given by Takahashi (2001) and 
Takaliashi & Ueliara (2001). By contrast, in an equiparti- 
tion distribution function 1/3 of the binding energy is sta- 
tistically distributed as internal excitation (Black & van 
Dishoeck 1987). Both produce an ortho-para-ratio (OPR) 
that is nearly 3. 

H2 is formed by associative detachment from H~ in a 
dust-free gas. This process is important in the clouds with 
lower dust content considered below. This is also an exother- 
mic process and we use the state-specific formation distri- 
bution function given by Launay et al. (1991). 

H2 is destroyed mainly via the 8olomon process, in 
which the H2 molecule is irradiated by far UV (912A< A < 
I200A) radiation and is excited to higher electronic states. 



Apjjroximately 10 per cent of these electronically excited 
H2 decay to the ground state continuum and are dissoci- 
ated. The other 90 per cent populate the higher vibrota- 
tional levels of the ground electronic state. These cascade to 
lower vibrotational levels producing infrared emission lines, 
referred to as 8olomon pumping. Formation pumping on 
grains is only 10 percent as effective as Solomon pumping 
when 8olomon destruction is balanced by grain formation. 
Thus, the H2 populations are non-thermal if the electronic 
lines are optically thin and the 8olomon process is domi- 
nant (hereafter referred to as the optically thin case), while 
the H2 level populations may be in LTE at the local gas 
kinetic temperature if the electronic lines are optically thick 
and the Solomon process is slow (hereafter referred to as the 
optically thick case). 

Radiative decays between ortho and para states are not 
possible because of the different nuclear spin. However, ex- 
change collisions with H, H+, and interactions on grain 
surfaces (below a certain critical temperature) can cause an 
ortho-para conversion. The column density ratio of J=l and 
J=0 traces the kinetic temperature in a coUisionally domi- 
nated gas but may fail to do so in a Solomon-process domi- 
nated region (AbgraU et al. 1992; Sternberg & Neufeld 1999; 
Roy, Chengalur & Srianand (2005)). 

The formation rate on dust has the largest uncertainty 
among the many processes considered in our calculations. 
There are significant variations in this rate even in the case 
of the Galactic ISM (Browning et al. 2003). There are also 
substantial differences between collisional rates of H2 com- 
puted by different groups at low temperatures (S05). These 
uncertainties should be kept in mind while comparing our 
results with observations. 



2.3 H I spin temperature: 

It is commonly assumed that the N(H l) spin temperature, 
Ts, is equal to the gas kinetic temperature. The optical 
depth of the 21 cm transition is proportional to N(H°) / Ts, 
so is sensitive to the inverse of the gas kinetic temperature. 
The mean value of T., we report here is Tk with weighting 
by N(H°) / T. A separate paper (Shaw et al. 2005) discusses 
our treatment of Ts, and relationships between T, Ts, and 
H2 temperature indicators, in detail. In DLAs, T, is usually 
estimated using the integrated optical depth Tv in the 21 cm 
absorption line and N(H i) measured from Lya using 



Ts = 



N{H i)U 



1.823 X 101St„ ■ 

Here, /c is the fraction of the background radio source cov- 
ered by the absorbing gas. Thus, low 21 cm optical depths 
could either mean high Ts or low /c. Even for /c = 1, the 
derived temperature will be high if the mean N(H l) cov- 
ering the extended radio source is lower than the measured 
one along the line-of-sight toward the optical point source 
(Wolfe et al. 2003a). Tims, observations will constrain either 
the physical conditions or the projected H I surface density 
distribution of the absorbing gas. 

2.4 Fine-structure level population: 

The ionization potential of C'^ is close to the energy of the 
photons responsible for the H2 electronic band transitions. 
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So, C I lines may be sensitive to the conditions in the H I — 
H2 transition region. The fine-structure level populations of 
C I are sensitive to the gas pressure and the IR radiation 
field. Thus, populations of the excited fine-structure levels 
of C I allow an independent probe of quantities that control 
the physical conditions in the H I — H2 transition region 
and the temperature of the cosmic microwave background 
(CMBR) (Srianand et al. 2000; Silva & Viegus 2002). The 
column densities of excited levels within the ground term of 
C I, C II , Si II, and O l are all calculated as part of the gas 
cooling function. All excitation and deexcitation processes, 
collisions, line trapping, destruction by background opaci- 
ties, and pumping by the external continuum, are included. 
At high z the IR pumping is predominantly by CMBR pump- 
ing, although the diffuse IR radiation from grains also con- 
tributes. 

2.5 Cloud geometry: 

We envision the region where the absorption lines form as 
a layer or cloud of ISM gas exposed to several radiation 
fields. In keeping with much of the PDR literature, we as- 
sume a plane parallel geometry (Tielens & HoUenbach 1985; 
Draine & Bertoldi 1996; Wolfire et al. 1995, 2003; Wolfe et 
al. 2003a,b). We further assume that the gas has constant 
density for simplicity. 

In the absence of ongoing star formation, the mcta- 
galactic background UV radiation field, dominated by QSOs 
(Haardt & Madau 1996), will determine the physical condi- 
tions and abundance of H2. If there is ongoing star formation 
then a locally-produced stellar radiation field will also con- 
tribute. OB stars are very short lived, and so do not move 
far from their birthplace before dying. So, newly formed OB 
stars will be close to their parent molecular cloud through- 
out most of their lives. This geometry is assumed in the PDR 
references cited above. 

Our main goal is to understand the physical conditions 
in the components with H2 and C I. Only the total H I 
column density is measurable in DLAs and the H I column 
density in the H2 component is generally unknown. So, we 
consider clouds with three values of 7V(H l); 10^'^, 10^", and 
10^^ cm~^. Wo assume the gas metallicity to be 0.1 Zq and 
vary the dust-to-metal ratio in the range 0.001 to 0.1 (this 
corresponds to a range in k (as defined in Section 1) of 10~* 
to 10"^) of the galactic ISM. 

We consider three ionizing continua; the meta-galactic 
radiation field at « = 2, the direct radiation field from an O 
star, and an O star continuum that has been attenuated by 
intervening absorption. The first mimics the case in which 
there is no in situ star formation. The second is observed 
in galactic star-forming regions - the OB stars are close to 
the molecular cloud and an H II region lies between them. 
This will be called the stellar case below. The third would be 
similar to a diffuse ISM exposed to the galactic background 
starhght, and will be called the diffuse case from now on. 

Following the general practice in the PDR literature, we 
define the intensity of the incident UV radiation field using 
a dimensionless constant x (as defined by Draine & Bertoldi 
1996), 

= -^912^ 1 (2) 

^ 1.22 X 107 ^ ' 



Here, Ama is the energy density (ergs cm ^) of photons and x 
= 1 for the Galactic UV field defined by Habing (1968). Thus 
X provides the UV field strength in the units of Galactic 
mean UV field. 

We use the observed metallicity, depletion, H2 abun- 
dance, and fine-structure excitations of C I, C 11, and 
N{C i) I N{C 11) to constrain either the particle density or the 
intensity of the radiation field. The 21 cm spin-temperature 
and the level populations of H2 are used for consistency 
checks. 



2.5.1 Ionization and thermal structure : 

In this sub-section we demonstrate the need for a composite 
self-consistent simulation of the gas in order to deduce the 
correct physical conditions using a pedagogical example. In 
Fig. 2.5 we show the ionization and thermal structure of a 
cloud irradiated by stellar and diffuse continua. 

Panel (a) plots the temperature of graphite and silicate 
grains (for the range of sizes considered in our calculations) 
as a function of depth from the illuminated side. All the 
calculations we present in this work use the self-consistently 
estimated grain temperature for a range of grain sizes which 
are important for different processes such as photoelectric 
heating and formation of H2 on the grain surfaces. 

The kinetic {Tk) and the electron density (rie) arc plot- 
ted in panel (b). In this work we present the H I 
weighted harmonic mean kinetic temperature as 
spin temperature (Ts). Detail investigations of relation- 
ship between Ts and Tk under different physical conditions 
is described in Shaw ct al.(2005). 

Panel (c) plots the densities of H2, C°, C+, and as 
a function of cloud depth. The ratio of carbon fine-structure 
levels are shown in panel (d). The electron temperature is 
high ('^10'' K) and the electron density is nearly equal to 
the density at the illuminated side of the gas for the 
stellar continua (panel b, solid line). A hydrogen ionization 
front occurs at a depth of 2.5x10^® cm, where T and Ue 
fall. Across the PDR, T ranges between 300 - 800K and 
electrons are mainly donated by C+ . The short-dashed lines 
show the results for the diffuse case. There is no H 11 region, 
and so the entire cloud is a PDR. The physical conditions 
are nearly constant across the cloud, which does not have 
enough grain opacity to attenuate the incident continuum 
significantly. 

The behavior of C° in the case of the stellar case is as 
follows. In photoionization equilibrium, 71(0°) is oc rie n(C"'") 
Qrec oc Tie n{C^) T~^'^ whcrc Qfrec is the recomblnatiou co- 
efficient. Here decreases by three dex across the ionization 
front, however the electron temperature changes by less than 
two dex. This leads to two orders of magnitude decrease 
in n(C°) across the ionization front (see upper panel). As 
the ionization potential of C" is very close to the energy of 
photons that are responsible for the excitations of electronic 
states in H2, one expects both C i and H2 to originate from 
the same part of the cloud. This happens for the diffuse case. 
But in the case of the stellar case, a considerable fraction of 
C I originates in warmer gas that docs not possess H2. 

The predicted ratio of fine-structure level populations 
depends on the nature of the radiation field. The ratios are 
constant across the cloud in the case of the diffuse case. 
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Figure 1. Pedagogical example: Various physical quantities are shown as a function of depth in a cloud irradiated by the stellar and 
diffuse continua with log Af(H i) = 20.7. The metallicity is 0.1 Zq, n}j= 50 cm~^, and log k is —1.39 (corresponds to a dust to metal 
ratio 0.4 times that seen in the Galactic ISM). In these panels the short-dashed lines are for the diffuse case. In panel (a) thick and 
thin curves are for Silicate and Graphite grains respectively. The short-dashed curves are for models with diffuse continua. The assumed 
radiation field is more or less 4 times that of Galactic mean UV field. In panel (b) the thick curves represent electron density. Panels (c) 
and (d) show the ionization and fine-structure excitations as a function of depth in the cloud. 



However, they strongly depend on the radiation field for the 
stellar case. 



UV field (i.e x = 1-44 x 10"^ for the Bgr at z = 2). Some 
of the results from our calculations are presented in Figs. 2 
and 3. 



3 IONIZATION BY THE META-GALACTIC 
UV BACKGROUND: 

First we consider the case in which the only available rar 
diation field is the meta-galactic UV background. We use 
the QSO dominated meta-galactic UV radiation field (Bgr) 
computed by Haardt & Madau (1996) and the cosmic mi- 
crowave background radiation (CMBR) at 2 = 2 (assumed 
to be a black body with T = 8.1 K). The UV flux density 
in the Bgr at energies below 1 Ryd at « = 2 is roughly two 
orders of magnitude lower than the current mean Galactic 
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3.1 Gas pressure: 

Panel (a) of Fig. 2 plots the moan pressure of the H I gas 
as a function of the total hydrogen density {nn) for three 
different values of Ar(H i). Although ionization and thermal 
gradients exist between the H 11 and H I regions, the H I 
gas is fairly isothermal. Thus, wo use the neutral hydro- 
gen weighted mean temperature to estimate the pressure. 
The continuous, short-dashed and long-dashed curves in all 
these panels are for dust to metallicity 0.1, 0.01, and 0.001 
(i.e K = 0.01, 0.001 and 0.0001) respectively . The regions of 
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Figure 2. The results for various constant density clouds ionized by the meta-galactic UV background given by Haardt & Madau (1996). 
The results with continuous, short-dashed and long-dashed are for k = 0.01, 0.001 and 0.0001 respectively. Panel (a) plots the mean 
gas pressure as a function of hydrogen density. The density-ranges of the cold and warm neutral medium are marked in this panel. 
The numbers near each line give the assumed log N(il l). This panel is useful to identify the warm and cold components of the stable 
two-phase medium under pressure equilibrium. The shaded histogram gives the observed distribution in the systems with H2 detections. 
The non-shaded histogram in panel (d) gives the results for the systems without H2. In panel (f) the non-shaded histogram provides the 
distribution of upper limits. The observational data used are mainly from Ledoux et al. (2003) and Srianand et al. (2005). The horizontal 
short-dashed lines in panels (b), (e) and (f) are typical detection limits obtained in echelle spectra. The horizontal short-dashed line in 
panel (c) gives the expected value of the ratio 7V(C l*)/Af(C l) when CMBR at z = 2 is the only source of excitation. 



thermal stability occur for d{logP) / d{logn)> 0. The warm 
neutral medium (WNM) and cold neutral medium (CNM) 
of the two-phase medium are shown for reference. As an 
example, gas with 0.3 ^ nH(cm^'^) ^ 1 is thermally unsta- 
ble for N(H l)~ 10^^ cm"^ The gas will be in the stable 
WNM phase for 0.03^ ni/(cm~^) ^ 0.3 and in the stable 
CNM phase for l^C nH(cm~^) ^ 30. For reference. Fig. 6 
of Wolfire et al. (1995) shows a similar phase diagram for 
various metallicities and dust content. The allowed mini- 
mum and maximum pressure in the two-phase medium are 
higher in our case than in the typical galactic ISM for a given 
column density, mainly because of the low metallicity and 



low dust-to-gas ratio (see also Petitjean et al. (1992), Lizst 
(2002), Wolfe et al. (2003a,b); Wolfe et al. (2004)). The main 
motivation for plotting the phase diagram from our calcu- 
lations is to have a rough idea of the nature of the gas at 
different densities and to compare our work with published 
models. It is worth keeping in mind the fact that ISM is 
more complex than different phases in pressure equilibrium 
with one another. For example magnetic field, if present in 
DLAs, can provide confinement even if there is no thermal 
pressure equilibrium between different phases. Thus, we do 
not make any serious attempt to model the DLAs as two- 
phase systems. 
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3.2 H2 abundance: 

In this section, wo compare the predicted and measured 
iV(H2) values (Ledoux et ah 2003) to determine the physical 
conditions in clouds both with and without observed H2. 

3.2.1 Systems without H2 detections: 

First we consider the cases where H2 is not detected. The 
predicted H2 column densities as a function of hydrogen 
density {iiii) arc shown in panel (b) in Fig. 2. The horizontal 
short-dashed line gives the typical detection limit achieved 
in H2 surveys (i.e., iV(H2)= lO" cm'^). The observed iV(H2) 
is distributed uniformly between 10^'' and 10^^ cm~^ (The 
histogram in the left hand side in panel (b)). 

From this figure it is clear that for a given A^(H l), the 
column density of H2 is independent of k when the density 
is low and the gas is mainly the WNM. This is mainly be- 
cause in the low-density, high-tomperaturc gas, H2 is formed 
predominantly through the process due to the low dust- 
to-gas ratio. It is also clear that the hydrogen density has to 
be higher than 0.1, 1.0, 30 cm'^ for iV(H i) = 10^\ 10^°, and 
10^® cm~^ respectively in order to detect H2. In the pres- 
ence of an additional local radiation field (perhaps generated 
due to in situ star formation) these critical densities will be 
larger. Thus, if the gas in DLAs is mainly a stable WNM 
in ionization equilibrium with the Bgr, then the equilibrium 
H2 column density will be below the detection limit. This in- 
ference is independent of k since the H~ process dominates 
the H2 formation at low densities. 

Fig. 2 suggests that if DLAs have a thermally stable 
CNM then the equilibrium abundance of H2 is high enough 
for the molecule to be easily detectable whenever k is greater 
than 0.0001 or dust to metal ratio grater than 0.001 (The 
long-dashed curves in panel (b)). The H2 formation time- 
scale, which can be long, does not affect this result. A typ- 
ical time-scale for forming H2 with molecular fraction /hj 
is ~/H2/2i?ri(H''). Here n{¥f) denotes the atomic hydrogen 
density. According to .Jura (1975) R~ 3 x 10"" cm+^ s~^ 
in the Galactic interstellar medium. Scaling this value by k 
we have, 

5.025 X 10-^ fn, 
t = =- yrs. (3) 

Assuming that all the hydrogen is H", we find a typical H2 
formation time-scale of ~ 5 x 10^ yrs with k = 0.0001 and 
fua = 10~^ for uh ~W cm~^. The age of the cloud has to be 
less than 10^ yrs for us not to detect H2 with Af(H2)^ lO" 
cm~^ and k = 0.0001. The hydrodynamical time-scales or 
pressure readjustment time-scales in the cold gas are usually 
larger than this value (Hirashita et al. 2003) due to the low 
sound speeds. Hence the typical age of the clouds is expected 
to be larger. Thus H2 should be detectable in a CNM with 
K more than 0.0001 in the absence of any additional local 
radiation field. 

3.2.2 Systems with H2 detections: 

Now we focus on the systems with detectable H2. Ledoux et 
al. (2003) have shown that these systems usually have a high 
metallicity and dust content, i.e; O.IZq and log k > —2. 
If the gas originates from a stable CNM, our calculations 



predict 7V(H2) ^ 10^^ cm^^ (Panel (b) of Fig. 2) for N{R l) 
> 10^° cm~^. The observed iV(H2) is always smaller. Inter- 
estingly, the observed Ar(H2) values with re ^ 0.01 are only 
reproduced in a very narrow density range, one that is usu- 
ally thermally unstable in the standard two-phase models. 
This signifies that, for a uniformly distributed range of den- 
sities, a cloud with a random choice of uh, with iV(H i) 
in the range of 10^° to 10^^ cm"^, and with re > 0.01, will 
have either no detectable H2 or N{'H2) > W^^ cm"^. The 
probability of detecting H2 in the observed column density 
range of 10^'' to lO^** cm ^ is very low. Thus, in order to un- 
derstand the relatively low column densities of H2 in DLAs 
with clear detections, we need either the presence of an ad- 
ditional radiation field or for the cloud to be too young to 
produce H2 . The existence of a local radiation field with in- 
tensity of the order or higher then the Galactic mean field 
has been suggested by various authors while discussing H2 
in individual DLAs (Black et al. 1987; Ge & Bechtold 1997; 
Srianand & Pctitjcan 1998; Pctitjcan et al. (2000;2002); Gc, 
Bechtold & Kulkarni 2001; Ledoux et al. 2002b, Levshakov 
et al. 2002; Reimers et al. 2003). 

3.3 Fine-structure excitations of C I and C ii: 

In this section, wo compare the column densities predicted 
for C I and C 11 fine-structure excitations with the obser- 
vations. Many of these results will not agree leading us to 
conclude that an additional source of radiation must exist 
in these systems. The time-scales to establish the ionization 
and populations within the fine-structure levels of neutral 
atoms are faster than the H2 formation time-scale, so this 
provides an indicator which should be in steady state. 

3.3.1 C I absorption: detectability and degree of ionization 

The N{C i)/N{C 11) ratio is a good tracer of the fiux of 
photons driving the Solomon process in the cold neutral gas 
where H2 forms since the ionization potential of C" overlaps 
with the H2 electronic bands. However, A'^(C 11) cannot be 
accurately determined in DLAs since the C llA1334 line is 
usually saturated. Unlike C 11, the column densities of Si II 
and S II (which usually trace the same region as C II (see 
Fig. 2.5)) are accurately measured using transitions with 
low oscillator strengths. Si and S arc usually not highly de- 
pleted in DLAs (but see Petitjean et al. 2002 for a unique 
counter example). i\r(Si 11) can be used as a proxy to esti- 
mate N{C 11) (Srianand et al. 2000) by assuming the solar 
abundance of Si/C. Now we compare the predicted N{C i) 
and A''(C i)/A^(Si 11) with the observations. We plot N(C l) 
and N{C i)/iV(Si 11) as a function of nn respectively in pan- 
els (e) and (f) of Fig. 2. The predicted values of N{C i) are 
much closer to the detection limit in a very low-density gas 
(riH^ 0.01 cm^'^). As previously noted, our calculations are 
performed with Z = O.IZ©, and the systems that do not 
show C I absorption tend to have metallicity in the range of 
O.OIZ0 ^ Z < O.IZq (Srianand et al. 2005). Thus, the ab- 
sence of C I is consistent with a cloud having nn^ 0.01 
cm^"^, for the metallicity and A''(H l) typically measured in 
these systems. Such a cloud will have log N{C i)/N{Si 11) 
< —1.5 (Panel (f) of Fig. 2). This is consistent with the 
measured upper limits, shown as the non-shaded histogram 
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in panel (f), in most of the systems without detected H2. 
But N{C l) is more than an order of magnitude larger than 
the typical detection limit for n^r^ 1 cm~^. Thus, C I should 
be detectable in a high-density gas, even for low metallicity, 
in the absence of any internal radiation field. The predicted 
values of N(C l)/N{Si 11) in the high-density gas is usually 
higher than the observed upper limits (see Liszt 2002 and 
Wolfe et al. 2003a). 

Next wc concentrate on systems with detectable C I 
absorption. As noted in section 1, apart from only one case 
{zabs= 2.139 toward Tol 1037-270), all the C I detections 
are from DLAs that also have H2. These systems have metal- 
licities higher than we assume. Based on the detection of H2, 
wc expect these systems to also have a higher density. Dense 
clouds {uh^ 0.1 cm~^) produce a higher value of iV(C i) 
than observed. The measured ratio of JV(C i)/N{Si 11) in 
these components is much less than our predictions for a 
high-density gas. Clearly higher radiation field is needed to 
produce N(C l) /N(Si 11) as measured in these systems. 

The fine-structure populations of C i and C 11 can also 
test the high density requirement for components with H2 
detections. This will be discussed in the next subsection. 



3.3.2 C 1 fine-structure excitation: 

Here we use the observed N{C i*)/N{C i) ratio to constrain 
hh in the components with H2. This ratio is regularly used 
to trace the pressure in a neutral gas (Jenkins & Tripp 2001; 
Srianand et al. 2005). Wc plot N{C i'')/N{C i) as a function 
of hydrogen density in panel (c) of Fig. 2. The dotted line 
gives the expected value of the ratio if CMBR pumping at 
« = 2 is the only source of C I fine-structure excitation. The 
observed values of A^(C l*)/N{C l) (the histogram in the left 
hand side) are much higher than this, suggesting that col- 
lisions and UV pumping also contribute to the excitation. 
Most of the observed ratios are consistent with the predic- 
tions for Hh ^ 10 cm~^ for the considered range of Af(H i) 
and K. As noted above, for such a high-density gas, our cal- 
culations predict A'^(H2) and N{C l)/N{Si 11) higher than 
observed. We show below that the presence of an additional 
radiation field can reduce both of these. 



3.3.3 Excitations of C 11 fine- structure level: 

Like C I, C 11* is always detected whenever H2 is present 
in DLAs. However, it is also seen in a considerable frac- 
tion of DLAs without C I and H2 (Wolfe et al. 2003a; Sri- 
anand et al. 2005). The observed column density of C 11* 
can directly give cooling rate (when the optical depth of 
[C ii]A158 line is negligible) which can be used to constrain 
the star-formation rate once the physical conditions in th gas 
is known (Wolfe et al. 2003a; 2003b; 2004). C 11 can originate 
from the CNM, WNM, and an ionized gas (i.e., H 11 regions). 
Collisions with atoms are important for exciting C I*, while 
electron collisions axe important for the excitation of C II*. 
Thus, C II* is expected to be detectable in systems with a 
warm and/or ionized gas even if the high-density CNM is 
absent (see Lchncr et al (2004) and Srianand et al. (2005)). 

C II* is invariably detected in all the systems with log 
A^(H i)> 21 and log — 0.032© irrespective of the absorp- 
tion redshift. The observed column density is in the range 




log Hydrogen density (cm ^) 

Figure 3. The calculated H I spin temperature is sliown as a 
function of density with the meta-galaetie UV background dom- 
inated by the QSOs at z ~ 2. The results are presented for three 
different column densities of iV(H i)(102i, lO^^, lO^o, lO^i cm-2) 
and three different values of k (continuous, short-dashed and long- 
dashed curves are for k = 0.01, 0.001, and 0.0001 respectively). 



12.7<log N{C ll*)<14.0 for systems at 1.5<2ai,s<2.5. The 
typical upper limit is N{C 11) < 10^^ cm~^ (Srianand et 
al. 2005). The calculated N{C ii*)/iV(Si 11) ratio is shown 
in panels (d) of Fig. 2. The shaded and non-shaded ob- 
served histograms in these panels are for the systems with 
and without H2 detection respectively. The observed range 
of N{C ii*)/iV(Si 11) in systems with H2 detections suggests 
they originate in a high-density gas. This is consistent with 
our conclusion based on H2 and A''(C l*)/A'^(C l). 

The observed N{C ll*)/N(Si 11) ratio tends to be smaller 
in systems without H2. As we have mentioned above, most 
of these systems have total N(ll l) higher than 10^^ cm~^. 
Thus, the systems that show C II* will originate in clouds 
with n/f^O.l cm~^. Thus, we will require a radiation field 
in excess of the Bgr to suppress C l and H2 in a high-density 
gas. 



3.4 H I spin temperature: 

The thermal state of H I gas can be probed with the 21 cm 
spin-temperature (Tg). The harmonic weighted mean tem- 
perature, a proxy for T^, is shown in panel (a) of Fig. 3. It is 
clear that if DLAs originate in a low-density WNM gas, then 
the spin temperature will be ~ 8000 K. Thus, systems with 
no H2, C I, C II*, and 21 cm absorption are consistent with 
a low-density WNM in radiative equilibrium with the Bgr. 
The predicted spin temperatures are usually less than 100 
K for clouds with N(H l)> 10^° cm'^ in the CNM. Thus, 
we expect all DLAs to show detectable 21 cm absorption 
if they originate in a high-density CNM gas which covers 
the background radio source. Unlike C i or H2, the presence 
of an additional radiation field with hv < 13.6 eV may not 
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Figure 4. The effects of micro-turbulence & cosmic ray ioniza- 
tion. The solid curves give the result without turbulence and cos- 
mic ray ionization when Z = 0.1 Zq and k = 0.01. The short- 
dashed curves show the effects of 3 km s~^ turbulence. The long- 
dashed curves shows the effect of cosmic ray ionization (with a 
cosmic ray ionization rate of hydrogen 2.5 xlO~^^ The re- 

sults are presented for two values of ^(H l)(10^'', 10^"*^ cm"'^). 



reduce the 21 cm optical depth because it can not ionize 
hydrogen. The absence of 21 cm absorption in most DLAs 
suggests that they originate in a low-density warm medium. 
Our calculations also suggest that Hi, H2, and C 11* should 
be found in systems with 21 cm absorption. The absence 
of C I and H2 in the few systems with 21 cm absorption 
is inconsistent with this prediction, again suggesting a local 
source of star light. 



cm ^. We would need 6 S> 3 km s ^ to produce a significant 
eflfect on Af(H2). 

[C l] 610/i is optically thick at high column densities and 
the fine-structure level populations are influenced by line 
trapping. Turbulence increases the photon mean free path 
and reduces this trapping. The effect is clearly seen for the 
ratio Af(C i*)/Ar(C i) when N(Q, i*)^ 3 x 10^^ cm'^ which 



occurs when uh 
10 

Af(H i) 



10 cm"^ (Panel (c) of Fig. 4) for Ar(H l) = 
. The effect is not seen significantly in the case of 
10^" cm~^ since large column densities of N(G i*) 
only occur for densities above 100 cm~^. As pointed out be- 
fore, our calculations produce higher N(Q, l) than observed. 
^'(C l) is invariably not saturated in most DLAs and A610/i 
line trapping should not control the N(Q, l*)/N{C l) ratio. 
Thus the inclusion of micro-turbulence can not rectify the 
problems of the Bgr models in simultaneously reproducing 
the H2 and C I observations. 

The column density of C 11* is also affected by the pres- 
ence of turbulent motions for N{C 11*) ^ 3 x 10^^ cm~^ due 
to optical depth in the [C 11] 158/i line. The observed A''(C 11*) 
is always higher than 3x10^^ cm~^ (Table 1 of Wolfe et al. 
2003a; Srianand et al. 2005) in DLAs with N{}i l)^ 10^^ 
cm~^. Thus, line trapping effects may be important in pro- 
ducing the observed excitation of the C 11 fine-structure 
level. 

Cosmic rays add heat to a highly ionized gas and pro- 
duce secondary ionizations in a neutral gas. produced by 
cosmic ray ionization can cause ortho-para conversion and 
thermalize this ratio (Flower et al. 1994). We consider a 
cosmic ray ionization rate equal to the Galactic background 
ionization rate (~ 2.5 x 10~^^ s~^; Wilhams et al. 1998). 
These results are presented with long-dashed lines in Fig 4. 
The enhancement in the gas temperature increases the col- 
umn densities of C I* and C 11* in the high-density gas. 
The pressure of the neutral gas increases due to cosmic ray 
heating as expected. We also see a decrease in N(H2) at a 
given riH mainly because the increase in the gas temperature 
reduces the efficiencies in forming H2. 

Background cosmic ray ionization does not produce 
drastic changes for the low-density gas where the Bgr dom- 
inates (Fig. 4). A much larger cosmic ray ionization rate is 
needed to have an desired effect. 



3.5 Effects of micro-turbulence and cosmic rays 

We have neglected cosmic rays and turbulent motions in the 
models considered until now. In this section we investigate 
whether the inclusion of these processes help bring the Bgr 
models into agreement with observations. 

The presence of micro-turbulence increases the mean 
free path for line photons and reduces line centre optical 
depth. We consider a turbulent velocity corresponding to a 
Doppler b parameter of 3 km s~^. This shows the greatest 
possible effect since the typical measured b parameters of H2 
components are usually ^ 3 km s~^ (Srianand et al. 2005; 
Ledoux et al. 2003; Petitjean et al. 2002; Srianand et al. 
2000). 

We consider two column densities, iV(H l) = 10^^ cm~^ 
and N(n l) = 10^" cm"^, along with k = 0.01, and Z = 
0.1 Zq. Some results are presented in Fig. 4. A''(H2) changes 
very little for A''(H l) = 10^^ cm~^. However, the molecular 
column density is slightly lower in the case of A^(H l) = 10^° 
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3.6 Summary: 

The main results for a cloud irradiated by the meta-galactic 
UV background radiation are: 

• The presence of the QSO dominated meta-galactic ra- 
diation field can maintain a H2 abundance lower than the 
detection threshold for 0.1 cm~^, irrespective of the 
dust content and N{}i l). The presence of any extra radi- 
ation field in addition to the meta-galactic radiation field, 
or a slower H2 grain formation rate, increases this critical 
density. Thus the absence of H2 in 85 per cent of DLAs is 
consistent with the low density models. 

• The detection of C I absorption is inevitable whenever 
our line-of-sight passes through the CNM. However, the den- 
sity range that produces the observed N{C i'')/N{C 1) ra- 
tio also produces N{C l) and iV(H2) higher than the ob- 
served values. An additional source of radiation with ener- 




Figure 5. The effect of density: The results of calculations of a cloud in the radiation field of a star with a surface temperature 
40,000 K (stellar case). The cloud has log A'^(H l) = 20.7, Z = 0.1 Zq. The long-dashed, continuous and short-dashed curves are for log 
ni]= 0.0, 1.0, and 1.7 respectively. The labels 1,2,3, and 4 in the short-dashed curves are for log(K) = —1.4, — 1.6, — 1.8, and —2.0 
respectively. We mark these numbers only for the short-dashed curves. The dependences on re has the same sense for other values of 
nu as well. The results are presented as a function of x (see Eq.2). In each panel, the observed distributions are given as histograms (see 
caption of Fig. 2 for details) . 



gies ^ 13.6 eV is needed to account for the low values of 
N{Y{2) and N{C l) in these systems. 

• Like H2, the absence of 21 cm absorption in most of the 
high- 2; DLAs can be naturally explained if DLAs originate 
mostly in the WNM. A low-density gas, corresponding to a 
WNM, has a very large spin temperate [Ta ^ 7000 K). If 
DLAs are dominated by such a gas then 21 cm absorption 
will not be detectable. 

• A high-density gas has Ts ^ 100 K. This produces 
strong 21 cm absorption along with high values of N{}i2) 
and N{C l). The fact that Af(H2) and N{C l) are not seen 
in the few systems that do show 21 cm absorption suggests 
that an additional radiation field is present in these systems 
as well. This is consistent with results of detail investiga- 
tions of individual systems available in the literature (see 
references given in Section 3.3.1). 

• The too-large Af(H2) and A''(C l) column densities pre- 



dicted at higher densities cannot be explained by micro- 
turbulent motions (up to 3 km s~^) or cosmic ray heating. 
However, the non-thermal motions affect the fine-structure 
level populations when the infrared lines become optically 
thick. The observed log N{C if) is ^ 13.5 for Af(H l)^ 10^^ 
cm~^. The effects of line trapping of [C 11] A158 becomes very 
important in these systems whenever the turbulent motions 
are small. 

In the following section we explore the possibility of using 
in situ star formation to prevent the formation of too-large 
iV(H2) in high-density systems. We wish to point-out that in- 
clusion of radiation from the Lyman Break Galaxies (LBGs) 
can increase the flux of Lyman Warner band photons in the 
UV background by upto a factor 10 (see Haardt & Madau 
(2001)). This will make the Bgr radiation roughly 5 times 
less than the Galactic mean fleld. However, Section 4 shows 
that the required radiation fleld is much higher than this 
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enhanced Bgr. In addition, this enhanced UV field will also 
produce bimodel distribution of N(H2) contrary to what has 
been observed. 



4 IONIZATION BY YOUNG STARS: 

We add a stellar radiation field using a 40,000 K Kurucz' 

model atmosphere in addition to the metagalactic radiation 
field discussed above. We consider direct stellar radiation 
and stellar continuum attenuated by N(H i) = 10^° cm~^. 
These two cases are denoted by "stellax" and "difi'use", as 
discussed in Section 2.5. We showed that the meta-galactic 
UV background is sufficient to suppress the formation of H2 
in low-density gas (i.e nn^ 0.1 cm""'). Thus, in this section 
we mainly concentrate on the high-density gas needed to 
account for the observed N{C i*)/N{C i) ratio. We consider 
clouds with three values of iV(H i) (10^°, 10^° '^, and 10^^ 
cm~^), three values of n/f(l, 10, and 50 cm~^), and four 
values of n (10^^ 10"^ ■^ 10"^■^ and 10"^ *). In all of these 
calculations we assume Z = 0.1 Zq, turbulent velocity 6=3 
km s~^. Cosmic-ray heating is not considered in the models 
described below. 

We present the results of model calculations as a func- 
tion of X (as defined in Eq. 2). The calculations for stellar 
case will have more high-energy photons than the diffuse one 
for the same value of x- Therefore, for a given Ar(H i), the 
gas will be hotter and more ionized at the illuminated side 
of the cloud for the stellar case (sec Fig. 2.5). We present 
the results of our calculations for the stellar and diffuse case 
for iV(H i) = 5 xlO^° cm"^ in Figs. 5 and 6 respectively. 
The effects of changing iV(H i) are shown in Figs. 7 and 8. 

4.1 H2 abundance: 

Panel (a) of Fig. 5 shows the predicted column density of H2 
as a function of x for the stellar case, log A^(H l)=20.7, and 
3 values of nn (long-dashed, continuous and short-dashed 
curves are for nH= 1, 10 and 50 cm~^ respectively). For each 
values of n//the results are presented for 4 different n. For 
a given nn and Xi a higher dust content k produces a lower 
iV(H2). Naively we would expect 7V(H2) to increase with in- 
creasing K due to an enhanced probability of H striking a 
grain, and increased shielding. However, the gas temperature 
increases for higher k due to additional grain photo-electric 
heating (see the next section). This reduces the II2 forma- 
tion rate as shown in Fig. 1 of Cazaux and Tielens (2002), 
so, ^(H2) decreases. 

The calculations reproduce the observed range of N{'H2) 
for K X sJlOO, riH ~ 10-50 crn"'\ and log N{H l) = 20.7 
for the stellar case. The range becomes 1^ x ^300 for the 
diffuse case (Panel (a) of Fig. 6). We find that with nH=l 
cm~^, II2 should be detectable when x ^ 3 and x ^ 10 for 
the stellar and diffuse case respectively. Clearly, for a mod- 
erate local radiation field, x — 1 ^ 10, II2 will be detectable 
for m^l cm"^ and iV(H l) ^ 5 x 10^° cm"^. Higher m is 
needed to produce detectable iV(H2) for a lower Ar(H i) (Pan- 
els (a) in Figs. 7 and 8). The range of x that is consistent 
with the observed range of H2 is summarised in Table 1 
for all the scenarios discussed in this work. Observations 
of atomic fine-structure lines will further narrow down this 
range. 



4.2 Spin temperature and 21 cm optical depth: 

Panels (b) of Figs. 5 and 6 show the predicted Ts for 
log N(H i) = 20.7, 3 values of (Long-dashed, continu- 
ous and short-dashed curves are for nH=l,10 and 50 cm~^ 
respectively) and 4 values of k. It is clear from all these 
panels that for a given nn (curves with similar line style), 
Ts increases with increasing k mainly due to photoelectric 
heating by dust grains (labels 1, 2, 3 and 4 on the short- 
dashed curve show models with k in the increasing order). 
We also notice that for a given n (say top-most curve with a 
given line-style) and x the models with higher jii/have lower 
Ts. This effect is very prominent in the stellar case. This 
is mainly because Ts is very large at the illuminated side 
of the gas in the stellar case. Panels (b) of Figs. 7 and 8 
show the results with nH= 50 cm^"* for 3 different values of 
N(H l) (long-dashed, continuous and short-dashed curves are 
respectively for log N(H l)= 20, 20.7 and 21) for the range of 
K. It is clear that for a given k (say top most curves for differ- 
ent line-styles), clouds with lower Af(H l) will have higher Ts 
and hence lower r„(21 cm). In the diffuse case (from Figs. 6 
and 8) we notice that for a given n^and k, Ts gradually 
increases with x and saturates to a constant for large values 
of X- The increase in Ts with lower value of x is the effect of 
photo-electric heating. At larger x the grains become highly 
charged and total heating rate will become independent of 
X (Bakes & Tielens, 1994; Weingartner & Draine (2001a)). 
Thus Ts becomes independent of x at large x- 

The range of predicted Ts in the range of x constrained 
by the H2 observations is summarised in Table. 1. This ta- 
ble also gives the expected 21 cm optical depth obtained 
using Eq. 1. It is clear from the table that for a high- 
density gas (i.e TiH^lO cm"^) and log N(H i)^ 20.7, 21 
cm absorption should be detectable with an optical depth 
of r„(21 cm)// ^ 0.27. Especially for high k, the 21 cm 
optical depth becomes as low as 0.05, even when na=50 
cm~^, for log N{U i) = 20. Thus, the absence of 21 cm ab- 
sorption in high-density (or low temperature) systems will 
either mean that the average N(H l) along the radio source 
is much less than N(H l) seen along the optical sight- line 
(Wolfe et al. 2003a) or that the actual N{ll l) in the high- 
density component is lower (Kanekar & Chengalur, 2003). 
Absence of H2 in the systems that show 21 cm absorption 
with low Ts will indicate a radiation field much higher (i.e 
X >> 1). Our calculations also suggest that it is possible to 
detect H2 with low r„(21cm) for moderate nnisee Models 
A & AE in Table. 1). For example, when log 7V(H i) = 20.7 
and uh— 1 cm~^, the expected spin temperature is high 
(2300-5600 K) and 0.05 ^ r„(21cm)// ^ 0.12. Thus, 21 cm 
absorption will either be weak or undetected in cases where 
H2 is detectable. Thus, the fine-structure excitation of C I 
or C 11 with systems with detectable 21 cm absorption will 
lead to a better understanding of physical conditions in the 
gas. This is detailed in the following sections. 



4.3 C I absorption: detectability and level of 
ionization: 

Panels (e) and (f) of Figs. 5 and 6 plot 7V(C l) and 
N{C l)/iV(Si 11) as a function of x for stellar and diffuse 
cases respectively with log N(R i)=20.7. It is clear from 
the figure that the predictions for different k (curves with 
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Figure 6. The results of calculations of a cloud in the radiation field of a star with a surface temperature 40,000 K and attenuated by 
N{B l) = 20 cm-2 (diffuse case). Rest are same as in Fig. 5 

Table 1. Results for clouds irradiated by starlight. A-E are the stellar case and AE-EE are the diffuse case 



Parameters 


A 


B 


c; 


D 


E 


AE 


BE 


CE 


DE 


EE 


log N(H 





20.7 


20.7 


20.7 


20.0 


21.0 


20.7 


20.7 


20.7 


20.0 


21.0 


(cm- 


^) 


1 


10 


50 


50 


50 


1 


10 


50. 


50 


50 


log N(H 


I(cxt)) 












20.0 


20.0 


20.0 


20.0 


20.0 


X 




^ 3 


1,30 


3,100 


0.5,35 


4,400 


^ 10 


1,100 


3,200 


0.3,22 


10,475 


Ts (K) 




3150.5600 


100,1000 


56,316 


104.1260 


48,275 


2290,5248 


90,900 


54,339 


40,560 


66.195 


T^{21 cr 




0.09,0.05 


2.74,0.27 


4.90,0.87 


0.53,0.04 


11.42,2.00 


0.12,0.05 


3.04,0.30 


5.07,0.81 


1.37,0.10 


0.83,0.28 


log N(C 


l)/N(Si 11) 
I*)/ N(C 
ii*)/N(Si II) 


-2.6, -2.0 


-2.6, -2.0 


-2.5, -1.8 


-1.8, -1.0 


-3.0, -2.5 


-3.8, -2.5 


-4.0, -2.04 


-3.7, -2.0 


-2.5, -1.0 


-3.7, -2.0 


log N(C 


-0.76,0.0 


-0.3,0.3 


-0.2,0.4 


-0.1, 0.44 


0.0, 0.40 


-0.76 


-0.6, -0.5 


-0.4, -0.1 


-0.4, -0.1 


-0.26, -0.44 


log N(C 


-l.S, -1.3 


-0.5, -1.0 


-0.2, -0.8 


-0.5, 0.0 


-0.9, -0.1 


1.9 


-1.3, -1.1 


-0.9, -0.6 


-1.5, -0.5 


-1.08, -0.66 


log N(C 


l(tot)) 
I*) 


12.0,13.4 


11.9,13.0 


11.8,12.8 


11.0,13.4 


11.8,13.0 


11.0,13.0 


11.0,13.0 


11.0,13.0 


11.2.13.0 


10.6,12.0 


log N(0 


11.2,11.3 


11.3,11.9 


11.5.12.2 


11.0,11.3 


10.5,11.7 


11.0,11.2 


10.0,11.5 


10.0,12.0 


10.0.11.5 


10.5,12.0 


log NCO 


I**) 


11.2,11,3 


11.3,11.9 


11.6,12.2 


11.0,11.3 


10.5,11.5 


11.0,11.3 


11.0,11.53 


110.0,11.5 


0.0,11.3 


10.0,11.7 


log N(Si 


II*) 


< 10.2 


10.5,11.0 


10.8,11.7 


r-^ 9.5 


9.6,10.2 


9.0.10.3 


9.9,10.0 


9.0,10.5 


9.0,10.6 


9.0,10.5 



same line-style) are very much identical. Both N{C l) and 
N{C l)/N{Si 11) are lower than those produced by the Bgr 
(Fig. 2) due to the presence of additional ionizing photons. 
A minor reduction in N{C l) and N{C l)/N{Si 11) for a given 
uh and k is noted for the diffuse case compared to the stellar 
case. This can be easily understood using Fig. 2.5. N(Si 11) is 
nearly identical for both the radiation fields. Whereas N(C l) 



is less in the diffuse case as and T are lower in these mod- 
els compared to stellar case (see panel (b) of Fig. 2.5). 

Here we concentrate on systems with H2 detections. 
The predicted 7V(C l) is well below the detection limit for 
log A'^(H l) — 20.7 and x J? 10 for the range of uh and the two 
stellar continua considered here. Thus, these models can ex- 
plain the weak or non-detection of C I absorption in some of 
DLAs that show strong H2 absorption {zabs~ 3.025 toward 
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Figure 7. Effects of column density: The results of calculations of a cloud in the radiation field of a star with surface temperature 
40,000 K (stellar case). The cloud has n^=50 cm~^, Z = 0.1 Zq. The long-dashed, continuous and short-dashed curves are for log Af(H l) 
= 20.0, 20.7 and 21.0 respectively. The labels 1,2,3, and 4 in the short-dashed curves are for log(K) = —1.4, — 1.6, — 1.8, and —2.0 
respectively. We mark these numbers only for the short-dashed curves. The dependences on k has the same sense for other values of 
Af(H l) as well. The histograms are as explained in Fig. 2. 



Q 0347-383 with log Af(H l)= 20.56; z^t,s= 2.595 toward 
Q 0405-443 with log Af(H l)= 20.90; and z^bs^ 2.811 to- 
ward Q 0528-250 with log Ni^ i)= 21.10). N(Q, i)/Af(Si ii) 
is lower than the measured ratio (in DLAs that show both 
H2 and C I absorption) for the range in x allowed by H2 
for both the diffuse and stellar case, and log Af(H l) = 20.7. 
Based on the trend seen in Figs. 5 and 6 we may need 
Hij^ 50 cm""^ and x ^ 10 to explain the observed range in 
iV(C l)/Af(Si 11) for H2 components that show detectable C I 
absorption. However, such a model will over produce Af(H2). 
This inconsistency can be solved by using a lower value of 
iV(H l) and a higher nH(see Panel (f) in Figs. 7 and 8). In 
these plots long-dashed, continuous and short-dashed curves 
are the results for n//=50 cm"^ with log N(H l) = 20.0, 20.7, 
and 21.0 respectively. The total Af(H l) measured for Zahs~ 
1.968 toward Q 0013-004 and Zabs= 2.087 toward 1444+014 
are consistent with log A^(H l)^20 in the H2 (and C l) com- 
ponents. The Zabs= 1.962 system toward Q 0551-366 that 



shows three H2 components has a total log A'^(H l)=20.5. 
The Zats^l.973 system toward Q 0013-004 has 15 weU de- 
tached C I components with a total N{}i l)=20.8. Clearly a 
low N(H l) is probable in the components that show H2 and 
C I absorption. A cloud with a moderate radiation field (i.e 
X ^ 10), log A'^(H l)=20.0, and n/f= 50 cm~^ reproduces the 
observed range of N{C l)/N{Si 11), 7V(H2) and N{C l). 



4.4 C I fine-structure excitation: 

In this section we consider the fine-structure excitation of 
C I. Panel (c) of Figs. 5 and 6 shows N{C i*)/N{C i) as 
a function of x for log N{R l)=20.7. This ratio should be 
independent of x if UV pumping is negligible because Ts 
depends mainly on density and is roughly independent of X- 
This happens for the diffuse case (Panel (c) of 6). Also the 
effect of K is clearly evident in this case. However, in the stel- 
lar case, we find that the predicted N{C i*)/N{C i) increases 
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Figure 8. The results of calculations of a cloud in the radiation field of a star with surface temperature 40,000 K attenuated by Af(H l) 
= 20 cm~^, our diffuse case. Rest arc same as in Fig. 7 



with increasing x- We also notice that for a given uh and 
X, the stellar case (Panel c in Fig. 7) with a lower Af(H l) 
produces a higher value of 7V(C l*)/7V(C l). However, the 
dependence on Af(H l) is very weak in clouds irradiated by 
the diffuse radiation field (see panel c in Fig. 8). This implies 
that the ratio 7V(C l*)/N{C l) will depend only on density 
for the diffuse case. However, the ratio 7V(C l*)/N{C l) will 
depend on the strength of the radiation field (also see panel 
(d) in Fig. 2.5) in the stellar case. 

Now we focus on systems that show detectable C I and 
H2. The predicted value of 7V(C l*)/A^(C l) is more sensi- 
tive to riH and weakly depends on A'(H l) and x for the 
diffuse case. The observed A''(C l*)/N(C l) is consistent with 
Wi^nn^ 100 cm"^. In the stellar case N{C i*)/N{C i) de- 
pends on Uh, N(H i) and X- Clouds with log 7V(H l) = 
20.7 reproduce the observed range in N{C l*)/N{C l) for 
l^n_f/SC50 cm~'^ and the values of x constrained by the H2 
observations. As noted before, however, these models fail to 
reproduce the observed A''(C l)/N{Si 11). Thus we require 
low 7V(H i) (~ 10^° cm-^), high nif(~ 50 cm"^), and low 



X (^ 10) components in order to be consistent with the 
observed Ai'(C i) and N{C i)/N{Si 11) ratio. 



4.5 C II fine-structure excitation: 

Here we discuss the predicted C 11* in detail. Panel (d) of 
Figs. 5 and 6 shows 7V(C ii*)/N (Si 11) for log N(}i i) = 20.7 
as a function of x- The shaded histogram is the observed 
distribution of the systems with H2 components and the 
non-shaded histogram represents those systems that do not 
show detectable H2 and C I. 

In the stellar case N{C ll*)/7V (Si 11) is higher for higher 
X (see panel (d) in Fig. 5), and for a given x the excitation is 
more for lower density (long-dashed, continuous and short- 
dashed lines are for nH= 1, 10, and 50 cm~"^ respectively). 
The ratio depends only weakly on k (different curves with 
different line- style) in the stellar case. All these trends are 
mainly because, for a fixed A''(H l), a considerable fraction of 
C II will originate from regions where hydrogen is ionized. 
The fraction of C 11 originating from a hot ionized gas is 
higher in the case of lower uh and so the ratio is higher. 
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Figure 9. The ratio of column densities of H2 in different rotational levels are sfiown as a function of total H2 column density for the 
stellar case. The points in the figures give the observed data (Ledoux et al. 2003). In all these calculations we assume the metallicity to be 
0.1 Zq, log N(H l)=20.7, and log k is varied between —2.0 and —1.4 (different curves with same line-style). The long-dashed, continuous 
and short-dashed curves are for ni]= 1, 10 and 50 cm~^ respectively. 



This happens for the lower A^(H l) case also (see panel (d) 
in Fig. 7 where long-dashed, continuous and short-dashed 
curves are for log N(H l) = 20., 20.7 and 21 respectively). 

Results for the diffuse case are summarised in panels (d) 
of Figs. 6,8. It is clear that for a given uh, N(H l), and k the 
ratio increases with increasing x when the x is small. How- 
ever, at larger x the ratio becomes independent of x- This 
is due to the grain heating saturation for highly charged 
grains at high x (see Weingartner & Draine, 2001a). This 
is also the reason for the lack of dependence of spin tem- 
perature on X (see previous section). Thus, at high x the 
ratios mainly depend on uh in the diffuse case. The mod- 
els presented by Liszt (2002) use the fitting function given 
by Bakes & Tielens (1994) and show a monotonic increase 
in N(C ll*)/N(C 11) with an increase in x- Weingartner & 
Draine (2001a) (see figure 15 in their paper) show that the 
fitting function given by Bakes & Tielens (1994) over pro- 
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duce the photo-electric heating at high x a-nd because of this 
there will be an increase in N(C ll*)/N(Si 11) with increasing 
X even at large values of As our treatment is very close 
to that of Weingartner & Draine (2001a) we clearly see the 
effect of saturation of photo-electric heating by dust grains 
in our models. 

First we will concentrate on the systems with H2 detec- 
tions. In the diffuse case the range in uh that is consistent 
with A'^(C l*)/A'^(C l) also reproduces the observed range in 
N{C ll*)/Af(S II). In the stellar case with log Af(H l) = 20.7 
the observed distribution of N{C ll*)/A''(Si 11) is consistent 
with Uh in the range of 10-50 cm~^. The models with lower 
A'^(H l) tend to produce higher A^(C ll*)/A^(Si 11) for a given 
Uh and k. A cloud with the low A''(H l) and high nn that 
are required to reproduce N(C l)/N(Si 11) will over produce 
N{C ii*)/iV(Si II). Thus the observed iV(C ii*)/iV(Si 11) seems 
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Figure 10. Same as Fig. 9 but diffuse ionizing radiation is used in these caicuiations. 



to favor the diffuse radiation field. Tfiis produces a model 
that reproduces all the other observations. 

Now we concentrate on systems without H2 detections 
but showing C 11* absorption. There are two possibilities for 
the absence of H2 in these systems: Either (i) the gas has 
lower density and so is partially ionized with a higher tem- 
perature or (ii) the gas is at a high density in a strong UV 
field. In the diffuse case the ratio N(C ll*)/N(Si 11) mea- 
sured in systems without H2 are consistent with uh in the 
range of 1-10 cm"'' (see panel (d) of Figs. 6). Srianand et al. 
(2005) pointed out that Al iii absorption seen in these sys- 
tems could be a useful indicator of the ionization of the gas. 
We will return to this issue while discussing the predicted 
N{Al lll)/A''(Al 11) ratio (see Section 4.7). 



4.6 Rotational excitation of H2: 

Here we focus on the H2 rotational excitation predicted in 
our calculations. 



4-6.1 The ortho-para ratio (OPR): 

The OPR indicates the kinetic temperature when the H2 
electronic bands are optically thick (i.e., log(Af(H2))^ 16; 
Tumlinson et al. 2002). Srianand et al. (2005) have shown 
that the OPR observed in DLAs are higher than those mea- 
sured in Galactic ISM, LMC, and SMC sight lines. Here, 
we probe the reason for this difference. Panel (a) of Figs. 9 
and 10 plots the OPR as a function of A''(H2) for the stel- 
lar and diffuse cases respectively with log 7V(H l) = 20.7. 
For optically thin H2, when the Solomon process dominates 
excitations of H2, (i.e log N{}i2)^ 16), the predicted OPR 
is close to 3 for clouds with log Af(H l)=20.7. However, the 
OPR is greater than 3 for log N{Y{2) in the range of 16 to 
18. At #(H2) > 10^* cm~2 the OPR traces the kinetic tem- 
perature of the gas. 

Sternberg & Neufeld (1999) show that the high value 
of the OPR seen for intermediate A''(H2) is because the elec- 
tronic absorption lines of ortho-H2 become self-shielded at 
smaller column densities than para-H2. Thus, ortho-H2 ex- 
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Figure 11. The ratio of column densities of H2 in different rotational levels are shown as a function of total H2 column density for 
diffuse case. The points in the figures give the observed data (Ledoux et al. 2003). In all these calculations, we assume the metallicity 
to be 0.1 Zq, nu= 50 cm~^, and log k is varied between —2.0 and —1.4. The labels 1, 2, 3 and 4 in long-dashed curves are for log(K) 
= —1.4, —1.6, —1.8 and —2.0 respectively. The long-dashed, continuous and short-dashed curves are for log Af(H l) = 20., 20.7, and 21 
respectively. 



ists while para-H2 is destroyed. Thus, we expect the OPR 
to be larger in the case of a higher radiation field when 
IG^log A''(H2)^18. For a given A''(H2) (in the intermediate 
range), the predicted OPR is higher for higher nn in mod- 
els with log N(H l)=20.7. We also notice that, for a given 
A''(H2) and nn, the models with lower A''(H l) produce a 
lower value of the OPR (panel (a) in Fig. 11). Thus, the 
observed OPR with 16^1og A''(H2)^18 will require a lower 
Af(H l) and higher nn- This is consistent with what we in- 
ferred based on the 7V(C l)/A^(Si 11) ratio. A detail observa- 
tion of an individual system confirms that the components 
with low OPR are consistent with a low value of A''(H l) 
(see Table. 1 of Srianand et al. 2005). As an example, Zabs = 
1.96822 toward Q 0013-004 has the lowest OPR value mea- 
sured in DLAs (0.64±0.09) and has log iV(H2) = 16.77 and 
log A''(H l)sC 19.43. We notice that the kinetic temperature 



is in the range of 40 — 560 K for the consistent models. 
This is slightly higher than the kinetic temperature range 
(60— 300K) derived using the OPR and assuming LTE as- 
sumptions (Srianand et al. 2005). We notice that the OPR 
does not track the kinetic temperature well in the intermedi- 
ate N(H2) range. A careful investigation of this is presented 
elsewhere (Shaw et al. 2005). All of our calculations are in 
qualitative agreement with the OPR seen for log A'^(H2)>18 
components. 

4.6.2 N(,J=4)/N(,J=2) and the radiation field: 

Panel (d) of Figs. 9 and 10 plots N(J=4)/N(J=2) as a func- 
tion of 7V(H2) for various log A^(H l) and k. The Solomon 
process controls these populations since the energy separa- 
tion between these energy levels is far too large and colli- 
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sional excitation is inefficient. This ratio indicates x when 
the H2 column density is low (Jura 1975). For a given iV(H2) 
with log 7V(H2)^ 16.0, the N(J=4)/N(J=2) ratio is larger 
for larger (i) uh, (ii)«: and (iii)A''(H l) (panel d in Fig. 11). 
Apart from the two systems in Ledoux et al. (2003), ab- 
sorption from the J=4 level of H2 is not detected. These 
two measurements and the upper limits for the optically 
thin systems are consistent with a radiation field as high as 
X = 30. There is very little difference between the diffuse 
and stellar continua since the excitation is mainly by elec- 
tronic fine absorption. In the optically thick cases, N{}i2) in 
the J=4 level is populated mainly by formation pumping. 
No clear trend is present since formation pumping depends 
on various quantities. 

4.6.3 N(J=2)/N(J=0) and N(J=3)/N(J=1): 

The N(J—2)/N(J—0) ratio is more sensitive to collisional 
excitation than the population ratios of higher rotational 
levels. The observed and predicted N(J=2)/N(J=0) and 
N(J=3)/N(J=1 ) are plotted as a function of N{\i2) in panels 
(b) and (c) of Figs. 9, 10, and 11 respectively. For the inter- 
mediate range of Af(H2) the models that reproduce the OPR 
also roughly reproduce these two ratios. However, they do 
not explain the observed distribution for log A^(H2)^ 18. It 
is important to note that the J=2 and J=5 levels are mainly 
populated by cascades from high J levels following formation 
and UV pumping. The grain formation distribution func- 
tion and grain surface interactions can affect the excitation 
of these high J levels. Although fitting the observed results 
may shed light on a gas with a different metallicity and dust 
composition than the Milky Way, such an exercise will divert 
us from our main theme and is left to future work. 



4.7 iV(AI lll)/Af(Al II): 

Here we focus on N{A\ lll)/A'^(Al 11) produced with a stellar 
radiation field on top of the Bgr(Fig. 12). Al 11 is mainly 
ionized by the high energy photons from the Bgr. In the 
diffuse case N{K\ lll)/Af(Al 11) ratio depends more on the 
density than on x- This is also the case for stellar case when 
X is small. However, the ratio increases with increase in x 
for large values of x (See thin curves in Fig. 12). Higher k 
produces higher 7V(A1 lll)/iV(Al 11), especially in the case of 
a high-density gas. The ratio Af(Al lll)/Af(Al 11) is higher for 
lower A^(H l) for a given k and nn. Our calculations predict 
log 7V(A1 lll)/Af(Al 11) to be less than —1.5 for the ranges of 
X, riH, and n that reproduce the observed properties of the 
H2 components. 

Srianand et al. (2005) have shown that most DLAs with 
log 7V(H l) > 21 show C 11* absorption even when H2 and 
C I are clearly absent. All these systems also show Al ill ab- 
sorption with log 7V(A1 m)/N[A\ 11) higher than-1.8 (Table 
6 of Srianand et al. 2005). Our model calculations produce 
log Af(Al lll)/iV(Al 11) higher than -1.6 for log 7V(H l)> 20.7 
when riH^ 10 cm~^. Clearly, the systems with only C 11* ab- 
sorption without H2 and C I absorption need nn^SlO cm~'^ 
in order to have a Af(Al lll)/N(Al 11) ratio at the detected 
level, since log A''(H l) ^21 in most of these components. This 
range in nnwill also explain the observed N(C ll*)/N(Si 11) 
in these systems (see Section. 4.5). This implies that C 11* 
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Figure 12. The ratio N{K\ lll)/Af(Al 11) as a function of x for 
clouds with log N(H l) = 20.7. The long-dashed, continuous and 
short-dashed curves are for n//= 1, 10, and 50 cm""^ respectively. 
As in the other figures, 1, 2, 3, and 4 marks the models with k 
— 1.4, -1.6,-1.8 and —2.0 respectively. The thin and thick curves 
represent the results for the stellar and diffuse case respectively. 

absorption originates in a region of lower density and higher 
ionization compared to the components that produce H2 and 
C I absorption. 

Lehner et al. (2004) show that a significant fraction of 
C II* absorption detected toward high latitude lines-of-sight 
in our galaxy originate from warm ionized medium (WIM). 
Using the profile coincidence between Al 11, Al ill, and pho- 
toionization models (computed using Cloudy) Wolfe et al. 
(2004) argue that C 11* in DLAs are unlikely to originate 
from WIM gas. They argue that a considerable fraction of 
Al III can be produced from the gas that is by and large neu- 
tral (much like the models considered here). For the Zaba = 
1.919 system toward Q 2206-19 that show C 11* without H2 
and C I absorption lines, Wolfe et al. (2004) derive a density 
of 1.6 cm"''. Srianand et al. (2005) discuss the profiles of 
absorption lines from different ionization states (including 
21 cm absorption line) in the case of Zaha= 1.944 system to- 
ward Q 11574-014 and conclude that a considerable fraction 
of C II* absorption originate from gas at lower densities (i.e 
riB— 1 cm""'). All these are consistent with our conclusion 
that density is lower (i.e 10 cm"'') in the systems that 
show C II* absorption without H2 and C I compared to the 
ones that show H2 absorption (i.e nH= 10-100 cm""'). 

4.8 Other fine-structure lines: 

We predict the fine-structure level populations of O I and 
Si II in addition to C I and C 11. The column densities of both 
O I* and O I** are in the range 2x10^^ - 10^^ cm"^ for the 
range of x suggested by the C I and H2 observations (see 
Tabid). The oscillator strengths of O I* and O I** lines are 
low (~ 4x lO"'^) and these lines are never detected in DLAs. 
Our calculations also predict N(Si 11*) < 2xl0'^^cm~^, con- 
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sistcnt with the fact that Si ll* absorption is not detected 
in DLAs. It is most likely that one may not detect these 
lines directly even if the absorbing gas has a higher density 
due to their weakness. However, it may bo possible to detect 
them by co-adding a large number of DLA spectra. Thus, 
one possible way to confirm the idea that most DLAs (with 
or without H2) originate in a high-density gas with star for- 
mation is to detect excited fine-structure lines of O I and Si 
II by co-adding many spectra as one does for metal lines in 
the Lya forest. 

4.9 Summary 

The main results for a high-density cloud in a stellar radia^ 
tion field are: 

• The observed properties of DLAs with H2 are consis- 
tently reproduced by models with a local radiation field 

in addition to the QSO dominated BGR. Most of the ob- 
servations of the H2 components (such as A*'(C l)/Af(Si 11), 
A^(C i*)/iV(C i) and A^(C n*)/iV(Si 11)) arc consistent with 
lower iV(H l) (i,e log TV (H l)~20 cm~^) and higher densities 
(10 ^ nif(cm"^) ^100). The median N{R i) in DLAs with 
H2 is ~ lO^'' * cm^'^, so in these systems only a fraction of 
the total A^(H l) originates in regions with II2. The typical 
kinetic temperature ranges between 40 and 560 K. 

• We reproduce the observed range of the OPR in DLAs. 
The systems that are optically thin in the H2 electronic 
bands have a lower OPR, suggesting log A'^(H l)~20, con- 
sistent with constraints from atomic species. The OPR > 3, 
seen in some of the components with intermediate H2 elec- 
tronic line optical depths, are produced by the different level 
of self-shielding in ortho and para H2. The absence of C I 
and H2 in J ^ 4 levels in the case of a few optically thick H2 
components are consistent with a higher x in these clouds. 

• Our predictions, the measurements, and the upper lim- 
its on the N(J=4-)/N(J=2) ratio in the optically thin H2 
components are consistent with a radiation field as high as 
X = 30. This is consistent with the limits on the radiation 
field from the atomic species. The absence of N(J=4) lines 
in the optically thick H2 components are consistent with a 
low rate of formation pumping in those systems. 

• H2 and C I are not detectable if the radiation field is 
much higher irrespective of the model parameters. However, 
such clouds will be easily detectable in 21 cm absorption 
with spin temperature in the range of 100 to 1000 K. These 
clouds will also show very strong C 11* absorption. However, 
the column density of C 11* will strongly depend on the 
amount of ionized gas along the line of sight. Also these 
systems will show very strong Al ill absorption. 



5 DISCUSSION AND CONCLUSIONS: 

5.1 Nature of the radiation field: 

Ledoux et al. (2003) show that detectable H2 absorption 
{N(H.2)^ 10" cm"^) is seen in 15-20 per cent of DLAs. We 
show that the observed properties of these systems arc in- 
consistent with a gas irradiated by the meta-galactic UV 
radiation field. Our calculations suggest that these systems 
originate from a high-density gas (^ 10 cm~^) irradiated by 
a moderate difi^use UV radiation field (1 to 30 times that of 



Galactic ISM) and indicate ongoing star formation in these 
systems. The mean radiation field is determined by both 
the SFR and radiative transport. As the mean dust optical 
depth in DLAs will be smaller than that of Galactic ISM, 
the typical SFR in DLAs with H2 absorption can not be 
much larger than that seen in our Galaxy. Even if such a 
moderate star formation exists in most DLAs, they will still 
contribute appreciably to the global star formation rate den- 
sity at higher rcdsliifts (sec Wolfe ct al. 2003a,b; Srianand 
et al. 2005; Hirashita & Ferrara 2005). 

5.2 Physical state of the H2 gas: 

Our calculations with a diffuse radiation field suggest high 
densities in the H2 gas (i.c 10 ^ nH(cm^'') ^ 100). The typ- 
ical temperature of the clouds that are consistent with the 
observations are in the range 40 to 560 K. Our calculations 
simultaneously explain H2 abundance and fine-structure ex- 
citations of atomic species without opting for enhanced H2 
formation rate on dust grains as required by analytical mod- 
els of Hirashita & Ferrara (2005). The inferred range in 
temperature and densities are consistent with the physi- 
cal conditions in the CNM gas. We show that if the cloud 
is irradiated by a difi^use interstellar UV background then 
N{C l*)/N(C l) can directly probe the density of the gas. 
For radiation field with x > 10 the ratio N(C ii*)/N(Si 11) 
will also trace the density of the gas as photo-heating satu- 
rates at higher values of x- However, if the cloud is close to 
the ionizing source N{C l*)/N{C l) depends also on x- The 
predicted values of N(C ii*)/N(Si 11) in the stellar case with 
low N(H i) (as required by other observations) are much 
higher than the observed values. Thus, our calculations re- 
quire that the H2 components in DLAs arc ionized by a 
diffuse radiation field. Most of the observations are consis- 
tently reproduced with N(H i) = 10^" cm~^. This suggests 
that only a fraction of the total measured N(H i) is present 
in the H2 components. 

5.3 DLAs without H2: 

The observations show that systems without detectable H2 
(i.e ~ 80 — 85 per cent of the DLAs) do not show C i ab- 
sorption and also have very small values of k. Roughly 50 
per cent of the DLAs show detectable C 11* absorption. Our 
calculations suggest that the absence of H2, C i, C 11* and 21 
cm absorption in a considerable fraction of DLAs could just 
be a consequence of a low-gas density in a moderate radiar 
tion field (irrespective of dust content of the gas). This more 
or less agrees with the Wolfe ct al. (2004)' conclusions that 
the systems with upper limits on C 11* absorption originate 
in a warm neutral medium (WNM)(also see Liszt 2002). 

The observed N(C n*)/N(Si 11) in systems that show 
detectable C 11* absorption without H2 and C i is consistent 
with nH> 0.1 cm~^ (see Section 3.3.1). The absence of C i 
and H2 in these systems can be explained as a consequence 
of higher radiation field. Wolfe et al. (2003a; 2003b) in the 
framework of stable two phase medium argue that most of 
the C II* absorption should originate from the CNM gas in 
order to have reasonable global star-formation rate density. 
We note nH<W cm~^ in systems that show C 11* without 
H2 and C I so that the observed N(C ii*)/N(Si n) as well 



(T^nnnnRAS MWRAsnnn nnn_nnn 



20 R. Srianand et al. 



as ionization state of Al can be consistently reproduced (see 
Section 4.7). This range is consistent with the one measured 
by Wolfe et al. (2004) for Zabs= 1-919 system toward 2206- 
10. Interestingly, the inferred density in these systems is less 
than that typically required to explain the property of H2 
detected components (i.e uh^ 10 cm^^). Thus, it appears 
that the systems that show only C 11* seem to originate from 
lower density gas compared to the ones that also show H2 
and C I absorption. 

Unlike H2 and C I, an additional radiation field (with 
\w ^ 13.6 eV) can not suppress 21 cm absorption in the 
high-density gas. Our calculations with rii/^ 1 cm~^ pre- 
dict a spin temperature in the range of 100 — 1000 K for 
a range of x, and A^(H l) typically seen in DLAs. Thus, 
21 cm absorption is definitely detectable. Over the redshift 
range that is similar to the range used by Ledoux et al. 
(2003) for H2 searches (1.9^Za6s*S 3.5), only two out of 8 
DLAs show detectable 21 cm absorption (Kanekar & Chen- 
galur 2003). The rest of these systems have a lower limit on 
the spin temperature in the range of 700-9000 K. Both of the 
systems with 21 cm absorption also show detectable C 11* ab- 
sorption. Detail investigation of one of these systems (zabs = 
1.944 system toward Q 1157+014) shows that C 11* originate 
not only from the gas responsible for 21 cm absorption but 
also from other components (Fig. 19 Srianand et al 2005). 
Clearly C 11* traces a wider range of physical conditions. 
There are few systems (c.g Za.bs= 3.387 toward Q 0201+11 
and Zabs= 3.063 toward Q 0336-01) that show detectable 
C II* absorption without 21 cm absorption. The derived up- 
per limits on spin temperatures in these systems will mean 
very low CNM fraction along the sight lines if the gas covers 
the background radio source completely (Kanekar & Clicn- 
galur, 2003). In the absence of VLBI observations, interpre- 
tation of these system will be very subjective (see Wolfe et 
al. 2003b for details). A careful analysis of Al iii/AI 11 and 
N(C ii*)/N(Si 11) in individual components is needed to get 
the contribution of ionized gas to the excitations of C 11*. 

Alternatively, one can use the fine-structure state pop- 
ulations of O I and Si 11 that trace C 11 very closely. Our cal- 
culations also compute expected fine-structure excitations of 
O I and Si 11. The expected column densities of O I*, O I** 
and Si 11* arc in the range 10^^ — 10^^ cm^^ for the range of 
parameters considered in our calculations. It may be possi- 
ble to detect these lines using pixel optical depth techniques 
that are used to detect metals in the diffuse low density 
IGM. Detection of such lines will put stringent constraints 
on density in these systems. 

5.4 Conclusions: 

In this article, we present calculations that self-consistently 
determine the gas ionization, level populations (atomic fine- 
structure levels and rotational levels of H2), grain physics, 
and chemistry. We show that for a low-density gas (nn^ 0.1 
cm~^) the meta-galactic UV background due to quasars is 
sufficient to maintain H2 column densities below the detec- 
tion limit (i.e iV(H2)^ lO^'' cm~^) irrespective of the metal- 
licity and dust content in the gas. Such a gas will have a 21 
cm spin temperature in excess of 7000 K and very low C I 
and C II* column densities for H I column densities typically 
observed in DLAs. 

Calculations with a high-density gas in the presence of a 



local radiation field reproduce most of the observations of H2 
components in DLAs. Thus our study clearly confirms the 
presence of CNM at least in 15-20% of the DLAs. We also 
show only fraction of total N(H l) is in the H2 components. 

Unlike the components with H2, interpretation of sys- 
tems that show only C II* without additional constraints is 
not clear. This is because presence of free electrons can be 
more efficient in populating the fine-structure level of C 11. 
This can lead to a high value of inferred tih if the electron 
contribution is neglected. Using Al ill absorption we show 
that a gas that produces C 11* in systems without H2 has 
lower density than the ones with H2 absorption. 
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